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1. INTRODUCTION 


One experimental difficulty »n measuring the low magnetic fields 
existing in the interplanetary medium Is that of having precise, 
absolute determination of the field vector. In the ideal case of a 
magnetically clean spacecraft and no zero-offsets in the instruments, 
the problem does not exist. In the case of a real S/C, there is 
always the possibility of a spurious field which includes the spacecraft 
residual field and/or a possible field from the sensors, due to both 

electronic drifts or changes of the magnetic properties of the sensor 

* • 

core (these latter effects may occur during the storage of the sensors 
prior to launching and/or in-flight). 

It is the purpose of this report to discuss some aspects of the 
problem and to test the reliability of a method originally devised by 
Hedgecock (1975). No mention will be given here about other methods 
which are described in the literature (see Davis and Smith (1968), 

Ness (1970), Ness et al. (1971), Rosenberg (1971), Belcher (1972), 
Neubauer and Schatten (1974). 

2. A BRIEF OUTLINE OF THE METHOD 

The basic postulate of the Hedgecock technique is that there should 
be no correlation between changes in measured field magnitude and changes 
in the measured inclination of the field with respect to anyone of three 
fixed cartesian component directions. By means of simple algebraic 
relationships between the field magnitude variation AB between successive 
measurements and the associated direction variation AP, one can show 
(Hedgecock, 1975) that the covariance function ^ between AB and Ap 
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for a given small spurious field correction c in the direction of a 
monoaxial sensor, is just proportional to the product c sin 0. 

This implies that a given field component, possibly affected by an 
unknown error c, can be simply corrected by searching for the value -c 
which makes the covariance function zero. In some experimental conditions 
it may not be so easy to find the zero crossing, due to the dependence of 
V^ c ^ upon the term sin 0: when o is small, it is possible that statistical 
fluctuations lead to more than one zero crossing, within the explored 
range of corrections. This ambiguity can be avoided by setting a minimum 
possible value for 0 reasonably larger than 0° (Hedgecock takes 30° and we 
do the same). 

There are a number of effects which may Influence the numerical 
value of the covariance function and hence the correction field, in a 
very complicated and practically unpredictable way. One source of error 
is the quantization of the experimental data. Another source -nay be 
identified when the set of data, or even part of it, used for the 
computation of is biased by a systematic correlation of AB and A©: 

it may happen for example that superimposed on a random variation there 
is a trend of B to increase during the period and of 0 to increase (or 
decrease). In some cases, when the technique is iteratively in succession 
for the three axes a bad determination of, for example the y axis, can 
mathematically lead to a significant second order correction which may 
well be physically meaningless. Some other sources are discussed later. 

In general, a few hundreds to a few thousands of data points prove 
sufficient for a good estimate of the corn ■ tion, so that the method 
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are smaller for shorter sampling times. Small variations of the spin 
rate do not significantly affect these conclusions. If data points are 
missing when averages are computed, a new source of noise occurs, whose 
effect is difficult to estimate If the pattern of missing data Is not 
known. The simplest way to reduce this source of noise Is to avoid 
using Intervals with too many mlsjing data points, 

4. APPLICATION OF THE MINIMUM CORRELATION TECHNIQUE 

4.1 A Simulation 

We took a field vector whose Intensity F and direction © and tf, where 
randomly variable around an average defined by certain values Fq, 

» we superimposed a spurious field defined by three orthogonal 
components c^, c , c^ on JF; and the minimum correlation technique was 
then used to determine the corrections. Table 2 shows some typical 
results. We used the same averaging Intervals considered above 
(At * 5 and 6 sec), together with 6 different sets of given corrections 
C l* C £» c 3 * ta bl e shows that when At * 6 the computed corrections 

are always within +0,27 of the specified corrections. When At = 5 sec, 
the discrepancy is small at the highest sanpling rate, but at lower 
sampling rates the discrepancies are large. Those discrepancies are 
drastically increased when the absolute values of c^, c ? , c^ are increased 
by a factor of two, which means that the method is very sensitive to the 
magnitude of the corrections. 

The occurrence of large computed corrections when At = 5 sec requires 
an explanation. The basic reason Is that using an averaging time which 
is not a common Integer multiple of the spin period and the sampling 
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period Introduces an artificial long-period modulation on the succession 


of *F' and 6F' superimposed on Its random part. This biases the 
x y 

covariance function, and as a consequence erroneous "corrections" are 
computed. Once a wrong value Is found, Its effect propagates to the 
following computations; In general, the process does not necessarily 
converge to any good value of the corrections. 


4.2 Application to (IMP-8) Data 



consecutive 15.36-sec-averages (treated as individual data). No 
significant difference was found between the three intervals, so the 
entire set of 72 zero-offset corrections was finally taken as a single 
set and statistically studied. 

The absolute value of the field correction on the X axis was found 


to be less than 0.27 in 44 cases and less than 0.37 in 58 cases, corre- 
sponding figures for the Y axis were <0.27 in 55 cases and <0.37 in 
66 cases. As regards to the Z axis, i.e. the spin axis, a slightly 
negative average field corrections of about -0.l7 was found. The 
results are sumnarized in Table 3 where average values and statistical 
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errors are shown, along with the linear best fit parameters (versus time). 

The results of Table 3 give us confidence that for real data the 
minimum correlation technique gives good estimates of the corrections 
Cj, c ,, c^, the only requirement being that data do not show a systematic 
correlation of AF and A0. 

A. 3 Application to Hellos 1 Data 

In this case the Individual data points that we used for the analysis 
are 5 sec averages. Early examination of the quick-look data showed some 
anomaly In ne In-flight offsets as compared to those measured In the 
pre-flight tests on the ground. This was at least partially due to some 
local spurious field from the spacecraft. 

Individual offset determinations have been made over each consecutive 
two-hour time interval. As an illustration of the results, we show in 
Fig» 2 a plot of the values of the field corrections versus time, 
together with the corresponding plot representing the flipper position. 

A striking ant icorrelat Ion shows up between the offset of the spin 
component and the flipper position. A simple interpretation is that there 
is a combination of both a S/C field and a sensor-associated drift. 
Individual contributions of the two sources must be determined by '-eans 
of some additional considerations. To this purpose, let us consider 
the situation illustrated in Fig. 3. By means of a flipper mechanism the 
sensors Dl and D3 are periodically rotated by 90° from the S/C spin 
direction to lie on the equatorial plane and vice-versa. The reference 
system X, Z is rigidly fixed on the S/C and rotates with it. 

Let us also assume that the spurious fields c^ and c^ can be 
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decomposed in two parts C., C, from the sensors package and S , S 

1 J x z 

from the S/C. 

In the position (a) the correction on z Is expected to be 

C - C, +■ *5 (3) 

J z 

As regards the equatorial component x, one can say that If the offset 
determined by standard techniques happens to be different from the 
nominal pre- flight value by an amount O' then 

D' - C l + S x (4) 

Similarly for the configuration (b) one gets 

C" - S z - Cj (5) 

D " - S x + C 3 (6) 

The above 4 equations then lead to the desired solutions 



which define completely both the S/C spurious field and the field 
perturbation to be attributed to the mechanical system which is rigidly 
rotating with the sensor itself. 

As an example we take from Fig. 2, C' * -0.37 and C" * +2.57. If 
the pre-flight determination of the sensors' offsets are compared with 
the actual flight determinations (see Villante and Mariani, 1977), we 
have for the time interval under consideration D' ■ 0,77 and D" ■ 2.17, 
and equations (7) lead to 
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2.87 


C 1 " * 2 * 17 S x " 

C- - -0.7 7 S - 0.47 

J i 

The fact that D' 4 0 and D" 4 0 means that some change occurred 

before launching or In flight, or both. One possibility Is that the 

change occurred In the sensors. Another possibility Is that at some 

stage of the post-test assembly of the S/C some magnetic component was 

used near ou- sensor. The only way to get ■ 0 would be C' ■ C" 

and D' ■ U' , In which case we would also have S ■ S ■ D'/2. On the 

other hand, the condlclon of no S/C field, S • S “0, tan be satisfied 

x z 

by a variety of situations, under the restrictions O' + C" ■ 0 and 
D" - C' "0, which simply Imply that Cj * -C" and ■ C'. In 
conclusion, it seems reasonable to believe that the corrections are due 
to the combination of a S/C field and a sensor package effect. 

In a later phase of the Helios l flight the corrections C' and 
C" changed significantly to become C' ■ 1.87 and C" * 27 which, 
combined with D' ■ 1.37 and D" ■ 1.07, leads to 


9 



C x • 0.05 7 


Cj - -0.27 

S x - 1.27 S z • 2.17 

Tills Is indicative of a temporal drift of both sources of spurious 
field, S/C and sensor package. 

It Is important to point out that our analysis Is limited to the 
Dl and D3 sensors. Whether or not a drift of the sensor 02 and a S/C 
field component parallel to It exist cannot be stated. The only 
Indication of a combined effect can be derived by the shift (If any) 
of the 02 offset from its nominal value, but there Is no way to decide 
which is the relative amount to be attributed to the twc possible sources. 
5. CONCLUS IONS 

The minimum correlation technique is generally appropriate for 
determination of the zero offset corrections of trlaxlal magnetometers, 
as investigated with S/C data taken at < 1 AU by IMP-8 and hellos 1-2. 

In general a number on the order of 1000 consecutive data points 
is sufficient for a good determination. 

The time separation between consecutive zero determinations can be as 
low as a few tens of minutes. However some care has to be exercised 
when individual measured points (rather than averages over appropriate 
time intervals) are taken, due to the possible unrecoverable effect of 
spurious field components perpendicular to the spin axis. The best way 
to avoid the difficulty is inat of using as individual points for the 
correlation averages over a time interval which is an integer multiple 
of the spin period and the sampling period. Use of simple algebraic 
relationships allows separation of the contributions due to the S/C 
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and the sensor system on those experiments with flippers. 
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TABLE 1 


TOE MAXIMJM VAU'ES OF AVERACE 6F x AND 5F y (IN /) 

WHEN At-5 SEC AND At-6 SEC, FOR DIFFERENT SAMPLING PERIODS T SA>fl>l 






SAMPLING 

TIME 


AVERAGING 

INTERVAL 

"CORRECTIONS" 

SPIN 

(SEC) 

0.1875 

(SEC) 

0.3750 

(SEC) 

0.75 

(SEC) 

1.50 

(SEC) 



1.00 

0.06 

0.14 

0.66 

0.67 



1.01 

0.10 

0.13 

0.64 

0.93 

At-5 SEC 


1.02 

0.13 

0.15 

0.60 

0.91 

|cj -2 \c 2 \ - 0 

1.00 

0.05 

0.10 

0.33 

0.67 




1.01 

0.07 

0.09 

0.45 

0.66 



1.02 

0.09 

0.11 

0.43 

0.93 



1.00 

0.00 

0.00 

0.00 

0.00 



1.01 

0.03 

0.04 

0.09 

0.13 

At *6 SEC 


1.02 

0.06 

0.07 

0.17 

0.26 

1 

I 

1.00 

0.00 

0.00 

0.00 

0.00 




1.01 

0.02 

0.02 

0.06 

0.09 



1.02 

0.04 

0.05 

0.12 

0.18 






































T\BLK 3 


Result* of the best linear fit ■ a^ + b of the temporal variation of 
the correction on the X, Y, Z direction as computed from IMP-8 data. 


Statistics on Algebraic Values 


(7 /day) 

bj (7) 

O (7) 

,00025 + .003 3 

-0.0183 + 0.0301 

0.25 

,00088 + .00264 

-0.0298 ± .0231 

0.19 

00677 ♦ .00325 

-0.109 + 0.028 

0.24 


Statistics on Absolute Values 


-.00079 + .00; '18 
-.00218 + .00178 
.00186 + 0.00198 












FIG. 2 



FIG. 3 



1. INTRODUCTION 


One experimental difficulty In measuring the low magnetic fields 
existing in the Interplanetary medium Is that of having precise, 
absolute determination of the field vector. In the Ideal case of a 
magnetically clean spacecraft and no zero-offsets In the Instruments, 
the problem does not exist. In the case of a real S/C, there is 
always the possibility of a spurious field which Includes the spacecraft 
residual field and/or a possible field from the sensors, due to both 
electronic drifts or changes of the magnetic properties of the sensor 
core (these latter effects may occur during the storage of the sensors 
prior to launching and/or In- flight). 

It Is the purpose of thl9 report to discuss some aspects of the 
problem and to t«st the reliability of a method originally devised by 
Hedgecock (1975). No mention will be given here about other methods 
which are described in the literature (see Davis and Smith (1968), 

Ness (1970), Ness et al. (1971), Rosenberg (1971), Belcher (1972), 
Neubauer and Schatten (1974). 



